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The zero temperature relaxational dynamics of a particle in a short range correlated random 
potential is addressed. We derive a set of "two-times" mean-field dynamical equations, accounting for 
a possible mean displacement of the particle when subject to an external force. We show first detailed 
results from the numerical integration of the above mentioned equations. We mainly pay attention to 
. the exponentially decreasing spatial correlations case, for which simple analytical arguments provide 

valuable results about the hessian (or the "instantaneous normal modes" structure) of the energy 
landscape, and we propose a geometrical description of the "mean- field aging". Our numerical 
results and further analytical arguments give access to the waiting-time dependence of the main 
characteristic time scales. 
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I. INTRODUCTION 



Understanding of the out-of-equilibrium dynamics of glassy systems, including spin glasses, structural glasses, 



supraconductor vortex glasses, etc, is a challenging problem. The need for exact, but non trivial results led to the 
introduction of ideal spin-glasses, like the celebrated Sherrington-Kirkpatrick model The spherical p-spin glass is 
such a model, where a closed set of equations exists for the time correlation functions of the dynamical variables, or 
soft spins 0(3). In the thermodynamical limit, each spin becomes coupled to an infinity of other spins. Each variable 
experiences gaussian fluctuations in the effective environment created by all the other spins. The dynamics then 
simplifies drastically, and reduces to a set of correlation functions which have to be determined self-consistcntly. The 
'"O ' result is the "dynamical mean-field solution" of the model. 

The mean-field dynamics of spin-glasses has revealed extremely rich, the most striking feature being the existence 
of a non-trivial aging relaxation regime at low temperature ||. For instance, the solution of the mean-field equations 
in this out-of-equilibrium regime demonstrate the existence of a generalised fluctuation-dissipation theorem (i.e. 
connecting correlation and response functions) whose validity seems now to extend to many realistic, non mean-field, 
models ||. 

Then, mean-field solutions are valuable for explaining the experimental aging of disordered magnetic systems . 
Finally, a sustained interest has followed the discovery of a deep formal analogy between the mode coupling description 
of structural glasses (supercooled liquids) and the mean field treatment of the spherical p-spin glass || . 

A crucial shortcoming of the mean-field description, however, is its inability to take into account properly thermally 
activated motion over energy barriers, leading to a sharp dynamical transition - divergence of an internal relaxational 
time scale - whereas the corresponding "finite dimensional" behaviour is only a strong but progressive slowing down 
of the dynamics. 

Despite of this last point, mean-field dynamics remains a major issue in the study of out-of-equilibrium statistical 
physics of disordered systems, and any approach providing a physical insight on its aging mechanism is of interest. A 
major step in that direction was made by J.Kurchan and L.Laloux M who investigated the zero temperature relaxation 
of systems including ferromagnets and spin glasses. The zero temperature limit makes it possible to consider the energy 
landscape, rather than a ill-defined "free-energy" landscape, without reducing the dynamics to anything trivial. 

In this work, we extend further their approach, and apply it to another system of interest: the mean-field dynamics 
of a particle in a short-range correlated random potential. 

The out-of-equilibrium, aging dynamics of this model has been first studied in H||, and thoroughly investigated 
in pc[] . Its glassy behaviour belongs to the same universality class than the spherical p-spin model. What makes this 
model interesting is its natural extension, when a finite and constant external force is applied to the particle. Then, 
it becomes a paradigm of "driven glassy system" , in which a non- linear response to the force as well as a significant 
violation of the fluctuation-dissipation theorem are expected, as shown by Horner |XT[ ] . 
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In this paper, we present the dynamical mean-field equations in the presence of a constant force, allowing for an 
arbitrary mean displacement along this one. These equations are then numerically solved, in the zero temperature 
limit, for an exponentially decreasing correlator, in the absence, and in the presence of a weak external force. The 
corresponding numerical results are presented, restricting ourselves to the linear response regime. Then, we start our 
geometrical analysis of the zero temperature relaxation by a simple random matrix calculation that we believe to 
describe satisfactorily the hessian of the exponentially correlated gaussian potential. Next, we perform an "instanta- 
neous normal mode" analysis of the relaxational motion. The key observable turns out to be the (intensive) energy 
difference between the energy £(t) at a given time t, and its asymptotic value lim^oo £(t). We subsequently analyse 
the waiting time dependence of two characteristic time scales tf, tf,, that we relate to £(t) — £(co). 

This work is preliminary to the study of the stationary driven situation in the presence of a finite force, which will 
be the subject of a forthcoming publication, and where the velocity-force characteristics, and the cross-over between 
linear and non-linear response will be exposed pjj] . 



II. THE OUT-OF-EQUILIBRIUM DYNAMICS IN THE MEAN-FIELD APPROXIMATION 



We introduce in this section the mean-field dynamics equations and discuss the low temperature aging solution in 
the absence of force. Let x(f) be the position of a particle, obeying a usual Langevin equation: 



x(i) = -W(x(t))+JP + C(i), 



(2.1) 



where are introduced the random potential V(x), the external force T, a white Langevin noise (,{t) corresponding 
to a temperature T, and a friction co effic ient equal to 1. Quantities x, V,.F, £ are A-dimensional vectors. Three 
sub-cases of the dynamics defined by (2.1) are of interest: 1/ the "isol ated " dynamics, without force: T = 0; 2/ the 
driven relaxational dynamics, which is the zero temperature limit of (2.1): J- =/= and T — > 0; 3/ the relaxational 
"isolated" dynamics: T = 0, T -> 0. 

The potential V(x) is a quenched disorder, chosen from a gaussian distribution. All the averages with respect to 
it will be denoted with an over-line • , while the average over the thermal noise (if any) £ will be denoted by the 
brackets ( • ). We suppose that the motion starts at t = and x(£=0) = 0. After averaging over the quenched disorder, 
this choice becomes equi vale nt to start with a random, "infinite temperature" distribution of initial positions. We 
expect that the process (2T) is self-averaging with respect to V(x) in the infinite dimensional limit. One introduces 
the correlator f{y) of the gaussian disorder, explicitly dependent on the dimension TV of the configuration space {x}. 



V(x) • V(x') — N ■ f 
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N 



V(x) = 0. 



(2.2) 



This form ensures a meaningful N — > oo limit, in which each coordinate Xj(i), or gradient component diV(x), remains 
of order 1, while the norms ||x(£)||, || VV|| scale like TV 1 / 2 . As a consequence, the external force must scale (ei being 
a unit vector) like: 



T = N 1/2 ■ F ■ ei. 



(2.3) 



One expects a displacement (x(£) ) = N 1 / 2 ■ u(t) ■ ei, and possibly a mean velocity (dx(i)/dt) = N 1 / 2 ■ v ■ e±. From 
now onwards, we arrange that ei coincides with the first coordinate axis i = 1. 

In the present paper, we restrict ourselves to the exponentially correlated potential: 



f(y) = exp(-y). 



(2.4) 



This is a special case of short range correlated random potential, characterised by lim^oo f(y) < oo. The average 
difference [V(x) — V(x.')] 2 is bounded when ||x — x'|| grows, and this ensures the existence of a normal diffusion regime 
at temperatures hi gh e nough. Another common choice is the power-law correlator: f(y) — 2/(7— l)-(l+j/)' 1-7 ^ 2 ; 7 > 
1 [p[-pl|. Choice ( |2.4| ) is also a particular case of f{y) = U 2 exp(— j//£ 2 ), with a pinning energy U p and correlation 
length £ set to 1, thanks to a simple rescaling, without loss of generality. 

The Langevin dynamics is handled with the help of a Martin-Siggia-Rose (MSR)-like functional integral, convenient 
for averaging over the gaussian disorder [ ^2|JT3]| . All the technical details corresponding to the saddle-point equations 
as N — » 00, are given in |l6| and the action is given in appendix The crucial point is that the limit N — > 00 
is taken first, before any other limit T — ► or t — > 00. As a result, we obtain a general effective quadratic action 
S[xj(t),ixj(t)], involving the original field Xj(t), and the MSR auxiliary field ixj(t). Three among the four following 
correlation functions appear explicitly in the action S[xj(t),iXj(t)}: 
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i(t) = N- 1 ' 3 { Xl (t)) ; (2.5) 



JV 



r(t,t') = N- 1 J2{xi(t)-ix j (t')h (2-6) 

3=1 
N 

b(t, = TV- 1 ]T ((^)-^(t')) 2 ); (2.7) 



J=2 
N 



d(t,t') = N- i j2((^(t)-^(t , w 

3=1 

= b{t,t') + [u(t)-u(t')} 2 . (2.8) 

These are the displacement u(t), the response function r(t,t'), and the correlation functions b(t,t') and d(t,t'). The 
Dyson equations for r, 6, d, u form a closed system of coupled integro-differential equations. For t > t' one has to solve: 

d t r(t,t') = 8(t-t') 

- [ ds 4/"(d(i, s)) r(i, a) [r(t, f) - r(«, *')]; (2-9) 
Jo 

d t b(t,if) = 2T- f ds4/'(d(M)) [r(t,s)-r(t',s)} 
Jo 

ds4f"(d(t,s))r(t,s) [b(t,s) + b(t,t')-b(s,t% (2.10) 

d t u(t) = F - [ ds 4/"(d(i, s)) r(t, s) ■ [u(t) - «(«)]. 
Jo 

(2.11) 

Equations (^9-2.11) are original ones, and allow for a non uniform displacement u(t). The aging, isolated, situation 
corresponds to the limit u(t) = F — 0, and d(t,t') = b(t,t') in the above system. The stationary limit, investigated by 
Horner |ll| amounts to write r(t, if) — R(t — t'), b(t, t') = B(\t — t'\), u(t) — v -t, and to reject the lower bound of the 
time integrals J dt ds to — oo. 

Three relevant observables: the energy £(t), the curvature M{t) and the pinning force F p (t) can be expressed with 
the help of these correlation functions. 



s^^n- 1 (y(x(t))>, 

rt 

ds2f'(d(t,s))r(t,s); (2.12) 

N 

3=1 

ft 

ds4f"(d(t,s))r(t,s); (2.13) 
F p (t)=N- 1 / 2 (-ckVW))), 

ds 4/"(d(t, s)) r(t, s) [u(t) - u(s)}. (2.14) 



The pinning force is such that (dx(t)/dt) = F p (t) + T with F p (t) = N 1 / 2 ■ F p (t) ■ e%. The pinning force F p and the 
driving force F have opposite signs. 

A proper study of the mean-field equilibrium phase diagram requires an extra quadratic confinement potential 
/xx 2 /2. This ensures the existence of a true thermal equilibrium in the high temperature phase, while correlation 
functions reach their asymptotic values exponentially fast. Then, a transition line Td(n), called dynamical tempera- 
ture, separates the high temperature ergodic phase, from a low temperature, aging and non ergodic phase p6|,p|,p[ . 

At high temperature, the system reaches a true stationary state, and the dynamics becomes time-translationally 
invariant (TTI), i.e. the 2-times correlation functions depend only on the difference t — t', while the 1-time expectation 
values are constant. In this stationary situation, it is convenient to introduce the TTI correlation functions B(t — t') = 
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limtjt'-xx) bit, t')\t-t' finite! R(t — t') — lim t .t> ^oo r(t, t')\t-t> finite- The fluctuation-dissipation theorem (FDT) holds 
and reads: 

dB(t)/dt = 2T- R(t). (2.15) 

As a consequence, the equal-time correlation functions coincide with their thermodynamical (canonical ensemble) 
counterparts. Taking the limit /i — > does not lead to any singular result |lC| ]. Provided the contribution from the 
harmonic potential has been subtracted off, the energy £(t) behaves smoothly as [i tends to 0. When [i exactly equals 0, 
the system cannot be at equilibrium, and instead, one expects a long time behaviour corresponding to a normal 
diffusion situation, with a finite diffusivity D = limt_ ) . o0 (x 2 (i) ) /(2Nt), a finite mobility rf 1 = limt— >oo u(t) / (Ft) , 
and the Einstein relation D — T ■ rf 1 . 

Kinzelbach and Horner described the dynamics in the stationary, high temperature phase [ ft6|| . They found that 
these correlation functions behave in the same way than those of the well known mode-coupling theories for supercooled 
liquids, as expected on general grounds |l7],[l8|]. The non-linearities of the self-consistent equations cause a dramatic 
slowing down of the dynamics as Td is approached from above, leading to a sharp transition at T = 

For instance, the function B(t), after a fast increase at short times t ~ 1, has a long plateau near a characteristic 
value B(t ~ tt) ~ q, before eventually reaching its asymptotic, long time regime B(t) = B(t/tb). Both tf and tb 
diverge like power laws of the difference \T — Td\ fl^l - 

The low temperature region however corresponds to an out-of-equilibrium situation. In the absence of external force, 
this is meant by the loss of both time-translational invariance (TTI) and fluctuation-dissipation theorem (FDT). The 
2-time correlation functions cannot be reduced any more to functions of the time differences t — tf, and there is a 
domain in the (t,t') plane, where the system ages [|8|— [l0| - 

The addition of a weak, constant external force leads to a somewhat different picture. As proposed by Horner, the 
system is expected to reach a stationary state (TTI), but the FDT remains definitively lost (TT). It turns out, however 
(cf next section) , that when the force is switched on at a time t = 0, there is a finite time interval during which 
the dynamics can be successfully described as a perturbation around the aging isolated (F = 0) regime, with a linear 
response approach. The extent of this linear response regime is inversely related to the magnitude of the force. 

The aging dynamics of the isolated particle has been exhaustively treated in |lC[] . The fluctuation-dissipation 
theorem is violated and must be replaced by: 

X{t,t')d v b{t,t') = r(t,t'). (2.16) 

In the time sector t — tf finite; t' — ► oo of the (t, t') plane, the behaviour is very similar to the one observed just above 
Td, and the value of X(t,t') is very close to its equilibrium value — 1/(2T). When the time separation t — t' ceases 
to be small relatively to a characteristic time t/(t') which has to be determined, X(t,t') departs from its equilibrium 
value, decreasing its magnitude \X\. 

The analytical study of the equations fl2.8| ) has only been possible in the asymptotic limit t,t' — ► oo, by dropping 
out sub-leading terms presumably of order 1/t, 1/t'. In this limit, the authors have shown (this is the crucial point) 
that it was possible to parametrise the dynamics with the help of the correlation function b(t,t') of the system itself. 
This implies that X(t,t') becomes a one variable function X[b(t,t')], and it turns out that all short range correlated 
models can be solved thanks to the ansatz ]Io| , f20|] : 

b(t,t') < q => X[b] = -1/2T; 

b{t,t')>q =>X[b]= X . (2-17) 

This extension of the FDT is called "quasi- fluctuation-dissipation theorem" (QFDT). In this paper, we rather use 
the function T(t,t') = —l/(2X(t,t')). In the aging regime, this effective temperature T = —l/(2x) is higher than 
the thermostat temperature T, and remains finite in the limit T — > 0. The physical meaning of these "two or many 
temperatures systems" is discussed in jgjj ]. 

In the case we are interested in, i.e. in the absence of confinement (fj, — * 0), for a correlator f{y) = exp(— y), \ and 
q are for any temperature T <Td, solutions of the system p0[: 




q e ~«/ 2 ; 



WW 



=9/2 



/2; 



and in the low temperature limit: 



T + T 2 /2 . 



T , ZT A 
4 ~r 16 



(2.18) 



(2.19) 
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In the same way, given a triplet t\ < ti < £3, in the time domain where (t\,t2,t 3 ) — > 00 and t\/t2, tijtj, finite, 6(^3,^2) 
is uniquely determined by the knowledge of b(ti,t\) and bjtp, t±). Again, the explicit dependence can be carried out 
exactly when the correlator is exponential. The result is Jig] : 

b(t 3 ,t 2 ) -q = b(t a ,h) -q- [b(t 2 , ti) - q]. (2.20) 

A well known shortcoming of this approach, is that any reference to the original times t,t' is definitively lost. 
The asymptotic solution cannot distinguish between b(t,t') and b(h(t), h(t')) where t 1— > h(t) can be any suitable 
reparametrization of the time variable. As a by-product, the previous analysis predicts only the more general form of 
the solution, in the aging regime t/t' ~ 1: 



b(t,t') = B 



hi 



fh(t) 



(2.21) 



\h[t>) 

For exponentially correlated potentials, the master function is known [ jlOj ] , and without loss of generality: 

b(t, t') = \n(h(t)) - \n(h(t')) + q. (2.22) 

In is made the conjecture h(t) — t 5 , compatible with the results found below. In what follows, we will refer to 
this solution as the time-reparametrization invariant (TRI) solution. 

At the b eginning of the aging regime, for t and t' such that (t — t')/t' is finite but small compared to 1, the scaling 
form ( 2.2l| ), reads: 



b{t,t') = B({t-t')-t^{t') + ...) +q. (2.23) 

Here, tb(t') is the characteristic time of the aging regime, defined by tb(t') = h(t')/h'(t'). This is the typical time 
needed by the particle for diffusing over a distance b(t,t') — q ~ 1. Non exponential correlators have a non analytic 
scaling function B(m) around m = and the r.h.s of ( [2.23 ) is singular in t — t' |1C 



The time-reparametrization invariant solution describes a situation where the time scales for the FDT regime 
(to = 1) and for the aging regime (tb(t')) are well separated (i.e. tb(t') ^ 1), which implicitly assumes t,t' — > 00. In 
order to go further, one needs to take into account the times derivatives d t , d t i neglected in the asymptotic regime of 
the TRI solution. It is enough, in principle, to fix up the reparametrization function h(t). Moreover, the TRI solution 
does not say how the parameter T goes from its FDT value (b(t, t') < q) to its QFDT value (b(t, t') > q). One defines 
for this purpose the new time scale tf(t), such that, for instance, T(t,t — tf(t)) takes a given value between T and 
— 1/(2%). We shall see below that tf(t) is much smaller than fa(t). 



III. RESULTS FROM THE NUMERICAL INTEGRATION OF MEAN-FIELD EQUATIONS 

The mean-field equations, with F and u(t) equal to zero were first numerically integrated by Franz and Mezard [^|,[|. 
The quadrature scheme is of order one in the time step h, but reveals itself surprisingly robust as h is increased up 
to value as large as 0.3. The authors of [|| report being able to reach t ~ 1000 at the best. Our investigations have 
shown that the quality of our solutions gets worst if h in increased above 0.2, and we present results up to t ~ 400. 



Fo r rea sons detailed in the next section, we have only considered the exponential correlator case (2^). We set T to 
in ( |2.8| ) and took the initial value C(0,0) = 0. The information coming from the numerics may be pigeonholed in 
three categories. 

1/ R esult s related to the TRI solution. First of all, we must check that the quasi-fluctuation dissipation 
relation ( 2.17 ) is true by plotting the integrated response versus the correlation function, on Figure (|). The observed 



value of x is close to 0.46, while the predicted value is 1/2. The TRI solution predicts also q ~ and lim^oo b(t, 0) 
00, in the absence of confinement. The measured asymptotic energy E (00) and mean curvature A'l(oo) are found to 
be in excellent agreement with the predicted values —2 and +4 respectively. 

2/ Beyond the TRI solution, without external force. This includes for instance the algebraic decay of 
the energy £(t) = — 2 + ci • t~ K . The exponent k is determined by plotting log(2 + £ (t)) versus log(i), and also by 
computing directly the logarithmic derivative, as shown on Figure (|). The exponent k lies between 0.66 and 0.67 
and our best estimate is c\ = 1.08. 

Also concerned are the characteristic times of the aging regime, and the precise nature of the cross-over from 
equilibrium to quasi-equilibrium fluctuation dissipation theorem. We are interested here in finding the characteristic 
time tf(t) as a function of t, defined by: 
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ds 2f'(b(t,s)) r{t,s) 



t-tf 



or alternatively: 



ds r(t, s) 



2/'(0) 



(3.1) 



(3.2) 



Equation ( jO] ) comes from the fact that the equilibrium X = — 1/2T and aging X = \ time sectors contribute for — 1 
each to the energy. The value tf which solves (3.1) separa tes the eq uilib rium regime (b(t,s) < b(t,t— tf)) fro m th e 
aging one (b(t, s) > b(t , t — tf)). The equivalence between (3T) and (3^2) is a straightforward consequence of ( 2.17 ) 
One generalises (3.2) in: 



ds r(t, s) = a, 



(3.3) 



For a < 1/2, t a must tend to a constant as t — > oo, while for a > 1/2, the asymptotic scaling (2.21) predicts: 



a-l/2 = - X B 



In 



h(t) 



h(t-t a ) 
X B(t a /t b ), 



(3.4) 



where terms (t a /tb) 2 have been neglected in the last expression, and ti, = h(t)/h!(t). If a is small enough, t a is simply 
proportional to Moreover, if hit) is indeed t s , then a — 1/2 = — x Bi((l — t a /tb)~ 5 ), and t a /tb is strictly constant. 
Our Figure (||) shows a plot of tf, i a =o.55, and £ a= o.45- The characteristic time scale tf tends asymptotically towards 
a power law C2t a , with C2 — 0.51 and a ~ 0.64 (according to our best estimate). 

The correlation function is found to grow logarithmically with t, and f(b(t, t')) = exp(—b(t, t')) behaves as a power 
law of t. Figure (||) presents exp(— b(t,0)) and exp(— b(t, t')) for a fixed t' . An algebraic decay t~ s of exp(— b(t, 0)) is 
likely, while exp(— b(t, t')) ha s not yet r eached its asymptotic regime, but could tend to the same t~ 5 behaviour. 

From the asymptotic form ( 2.2l| , |2~22] ) we note that exp(— b(t, t')) ~ h(t')/h(t); t, t' — + oo; t/t' finite, and our Figure 



is consistent with h(t) — t 6 , 5 ~ 1.10. Also shown is exp(— b(t,t — tf)) — t~ 7 , 7 ~ 0.42. Expanding exp(— b(t, tf)) 
as (t/tf)~ s , and using tf ~ t a , one finds a relation 7 = 6(1 — a) between exponents. The agreement between 
6(1 — a) = 0.39 and 7 ~ 0.42 is acceptable. 

3/ The linear displacement regime, in the presence of a driving force. A small force F is applied and the 
displacement u(t) monitored. The linear response implies that u(t) must be proportional to F, and it is indeed the case 
for time intervals not too large. Figure (||) presents u(t)/F for decreasing values of F. The curve F — 0.05 is virtually 
indistinguishable from the integrated response lZ(t) = f Q r(t,s)ds, and this shows that lim^^o UF(t)/F = TZ(t), 
i.e. the expected linear response behaviour. The other curves depart from the integrated response after a time tp 
decreasing with F. When starting from the isolated and aging situation F — 0, the linear response only holds during 
a finite time interval < t < tp. What happens later is the onset of a stationary state, with a well defined velocity 
v and a non- linear dependence in the force as advocated by Horner [0. A study of this regime is to be published 
soon p7[ . 

In the linear response regime, Figure (||) is compatible with: 



t(t) = F ■ (c 3 + c 4 • ln(*)), 



(3.5) 



C3 



0.71 and c 4 = 0.54. 



IV. THE GEOMETRICAL APPROACH 



In this section, we transpose to the particle in a random potential some of the ideas which have revealed fruitful 
when applied to the spherical p-spin model, namely the geometrical analysis of Kurchan and Laloux Q|. We expose 
first the main concepts of the method, and then propose a method for computing the limit value of the dynamical 
energy from generic properties of the potential V(x), working only for an exponential correlator f(y) = exp(— y). 

The p-spin model starts aging below a dynamical temperature Td , and encounters a thermodynamical glassy tran- 
sition at T s < Td p2[ . Detailed investigations have brought an appealing picture of the complex free-energy landscape 
of the spherical p-spin model, accounting for many features of its thermodynamics p3]-|25| and its dynamics |^,p^] . 
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The phase space of the p-spin model can be investigated with the help of a "Thouless- Anderson-Palmer" free-energy 
^(rrii) of the magnetization rrij, i = 1 . . . N. At low enough temperatures, the function $ develops many extrema m\ a \ 
the TAP solutions a. Those extrema which are minima, i.e. the second derivative matrix d 2 $> / 'dmidrrij is definite 
positive, are metastable states, as they are separated from each others by extensive free-energy barriers. A particular 
realization of the system, prepared in a given metastable state a remains for ever in this state in the thermodynamic 
limit. 

The stability of a metastable state is related to the lowest eigenvalue A m in of Spec(d 2 $/cmiiCmij), spectrum of the 

hessian matrix. \ m i n turns out to be a monotonically decreasing function of the free energy $(m| Q ^) of the state 
itself. This defines the free-energy <&d of the marginal states as X m in — for <f> = $d. Magnetizations such that 
$()7ij) > $<; represent regions of negative curvature which does not contribute to the thermodynamics but play a role 
in the dynamics [ fblf . 

The glassy dynamics of the p-spin model is observed when the stable metastable states $( Q ' < $<; are populated, 
or equivalently, when the canonical Boltzmann measure is split into its metastable components a. At a temperature 
lower than T^, thermal equilibration requires the system to explore all the relevant metastable states $ < <$>d- Such an 
equilibration is impossible as going from one metastable state to the next one requires to go over an infinite barrier. 
What happens instead to a system quenched from high temperature, to T < Td, is the onset of aging. The systems 
wander more and more slowly around the magnetization region $(mi) ~ $d, i.e. around the marginal states. 

The zero temperature relaxational dynamics is simpler because the free energy reduces to the hamiltonian H of 
the spins Si. At variance with the finite temperature case, the regions with negative curvature of H are now well 
defined. Taking the limit T — > in the mean-field equations does not lead to any singular behaviour. This somewhat 
counter- intuitive property is the consequence of sending N — > oo first, by keeping finite the times t and t 1 . The 
dynamics is a pure gradient descent, but remains non trivial. 

In order to perform a geometrical analysis of this relaxational dynamics, it is necessary to keep the dimension TV 
large but finite. Then, the relaxational process occurs till the particle falls into a true minimum of the Hamiltonian 
H(si), and, at T = 0, remains stuck there indefinitely. According to the description advocated in [Q, a system 
starting from a random configuration {si(0)} will explore regions with smaller and smaller gradient || Vi?(si)||, and a 
decreasing number of negative eigenvalues in the spectrum of Ti-ij, hessian of H(si). The typical time t J N for reaching 
regions where Tiij has only / negative eigenvalues, diverges as N is sent to oo by keeping / finite Q . As a consequence, 
in the N — > oo limit, the system is unable to reach within a finite time t a true minima, or even a saddle between two 
minima, and the difference £ — £d remains positive. 

To what extent does the above picture describe the particle in a random potential ? Much less is known about the 
properties of the metastable states, and there is no obvious equivalent of the free energy functional $(( x )) of the mean 
particle's position (x) ( - Q ' ) in the state a. Nonetheless, we expect that the basic mechanism of the dynamical transition 
remain the same as for the spherical p-spin model, i.e. a slow relaxation toward a region of marginal states, A mi „ ~ 0. 
When considering the zero temperature limit, the dynamics reduces to a gradient descent x(t) = — VV(x(i)). The 
metastable states now correspond to local minima of the potential V(x), and their stability will depend on the 
spectrum of the hessian Hij = Q?-F(x), where d{ means d/dxi and dfj = d 2 /dxidxj. 

The purpose of the geometrical approach, at zero temperature, is to relate the values provided by the more formal 
field-theoretical approach, to basic properties of the potential V(x). For instance, one must be able to compute the 
asymptotic values of the energy £ (t), curvature M(t), and the "plateau value" lim-r^o l/T of the correlation function. 
This is the first step, already outlined in [Q. One of the original contributions of this work concerns a second step, 
where we justify, with some geometrical arguments, many of the fine properties of the aging behaviour, beyond the 
time reparametrization invariant solution. 

A random matrix computation of the spectrum oiTL. A challenging problem in the study of the supercooled 
liquids dynamics, concerns the computation of the canonically averaged spectrum of the instantaneous normal modes 
( Spec(H) ) p5| . Here, TL is the dynamical matrix, hessian of the potential energy V of the interacting particles, 
function of the coordinates n. Spec(7i) is any representative characteristic function of the eigenvalues spectrum, e.g. 
the density of states. 

( Spec(H) ) = Z- 1 ■ [ dx ( e - fJV ■ Spec(H)) . (4.1) 

Z is the canonical partition function, and D is a bound domain eventually becoming infinite. In the mean field 
situation, as the energy V is exte nsiv e, the canonical average is dominated by a saddle point value V{0) of the 



potential V(x). The analogous of (4.1) becomes: 



Spec(H) ) = Wg 1 / dx tf(V(x) - V(fi)) Spec(H)(x), 



v 
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W = [ dx S(V(x) 
Jv 



(4.2) 



In the p-spin case, the spectrum of TL is self-averaging with x, i.e. the spectrum of Ti is a shifted semi-circle, 
by an amount controlled by V(f3) only. The dynamical energy is found to be the highest value of V{0) such that 
Spec(H) > (marginality). 

The averaged spectrum ( Spec(7i) ) defined by equation (4.2) seems to be the natural quantity to consider when 
looking at the zero-temperature relaxational dynamics of our mean field model. We have found that as far as expo- 
nentially correlated potentials are concerned , Spec(7i) is, at the leading order, a non-fluctuating quantity determined 
by Vq = V(P). More precisely, Spec(H) has got a semi-circular distribution of radius A, centred around V. 



A = 4VT(0); 
p _ 2/(0) Vq 
/(0) N ' 

Let us outline our demonstration. We consider first the (r independent) "annealed average". 



S(V(r)-V ) Spec(W)(r). 



(4.3) 
(4.4) 

(4.5) 



In order to compute (4.5), it is enough to enumerate the correlations of V (v) , dijV (v) , where r is an arbitrary point. 
All the dijV(r) are independent at the leading order A^ -1 , whereas the N+l remaining variables V(r),daV(r) are 
found to be correlated. One has: 



N-[d l3 V{v)f 



N ■ [d u V(rW 
N ■ daV(r) ■ djjVj r) 
V(r) ■ duV(r ) 
N- 1 • [V(r)] a 



4/"(0) 
12/"(0) 
4/"(0) 
2/'(0) 
/(0) 



-OiN- 1 ) 

-o(jv- x ) 

-OiN- 1 ) 

-o(jv- x ) 

-OiN- 1 ) 



(4.6) 



H. is split into a scalar part T>Sij and a fluctuating part 7i'. The elements of TC are independent and gaussian, and its 
eigenspectrum has, at the leading order, a semi-circular shape of radius Ay/ /"(0) centred around 0. If A~ — > oo and 
V(r)/N finite, then T> is constant, up to fluctuations of order N^ 1 / 2 (cf appendix [§]). 



V 



2f'(0) V(r) 
/(0) ' 



A^ 



0(N 



-l/2 ) 



(4.7) 



The resulting spectrum is the one ann oun ced in equations ( 4~3| , 4.4). 

In order to bridge the gap between ( |4.2| ) and (fl.5|), we consider now the two-points annealed average: 



S(V(r) - Vq) 5(V( ri ) - V x ) Spec(W)(r). 



(4.8) 



The analysis involves now correlations between V(r),dijV(r), V(ri),dijV(ri). One finds that for a generic correlator 
f(y), Spec(H)(r) depends on both V(r) and V(ri). However, if f(y) obeys f ■ f" — (/') 2 = 0, with f(y) = exp(— y) 
as a particular case, the dependence in V(ri) disappears, and the result (4.4) holds. 
Computing 



S( V(r)-V ) SiVir^-Vx)... 
x6(V(r n ) - V n )Spec(W)(r), 



(4.9) 



becomes very difficult as n > 3, and we were not able to find a close expression for Spec(H)(r) (V , Vi, . . . , V n ). 
However, if f • f" — (f 1 ) 2 = 0, again Spec(7i)(r) depends only on Vb, and (4.4) is valid. This shows that the spectrum 
of 7i is a local quantity, independent of the environment of the particle. 



Because Spec(H(r)) is a function of V(r) only, we conclude that the average fl4.2j ) is described by (f4.3| , [4.4! ) and that 
the self-averaging property of Spec(H)(r) and its linear dependence in V(r), which was true for then-spin model, is 
still true for exponential correlators. The appendix [b| gives further details on the computation of (4J5) and (4.S), 

Now, we suppose that the trajectory x(t) explores representative regions of the potential (i.e. non-excepti ona l 
points), for which the above mentioned results hold. The lowest eigenvalue —S of the hessian, defined by (4.4) 
becomes a time-dependent function: 
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S(t)=A-^-S(t), (4.10) 

leading to the energy dependent (through S) density of eigenvalues of Tiij . The number of eigenvalues between A — S(t) 
and A — S(t) + dA is p(X)d\ (time independent). 



p(A) = 2(ttA 2 )-V a (2A — A). (4.11) 



The marginality condition, by definition, is S = 0. Equation ( 4.1C ) yields the "geometrical energy", necessary for Ti 
to be marginal: 



and the curvature M geom : 



£ 2 vrmm u 12) 



Mgeom = I dAA p(A), 



= 4 v / 7 77 (0). (4.13) 

After a time t long enough, the particle evolves in a marginal region (S(t) ~ 0) of the potential V(x), with a 
small gradient ||VV(x)||. At low temperature, the potential may be developed up to the second order by means 
of local coordinates yc. V{y) — V{Q) + y • VF(0) + J2i=i n The plateau value "g" of the correlation 

function b(t,t') is thus given by assuming that each direction of curvature A^ is thermalized with / yf) ~ T/\ i: and 
q = 2N-^ l=hN {yf): 



Sy/f"(0) 2T 

T 



Qgeom = / PW\s=0, 



VTW) 

Let us compare now with the results from the dynamical mean-field theory, in the zero temperature limit |K 



(4.14) 



VJW) /(°) 

lim M(t) = 4 v / /"(0), (4.16) 

t—>oc 

(4.17) 



vfi 

Agreement holds for the curvature and q, whereas the geometrical and dynamical energy differ, unless /(0)/"(0) 
f'( Q) 2 . We cannot conclude about the relevance of the geometrical approach for a generic correlator, e.g. power law, 



as (4.4) probably does not hold. However, the exponentially correlated toy- model turns out to be a very favourable 
model, for which the geometrical approach gives reasonable results. The following of this paper aims at demonstrating 
that many features of the zero temperature dynamics of this model (exponents, aging, driving with a force) can be 
explained with the help of geometrical arguments. 



V. THE DISTRIBUTION OF THE GRADIENT'S COORDINATES 



In this section, we define an orthonormal frame "attached" to the particle. The procedure used is reminiscent 
from the definition of the instantaneous normal modes in the study of supercooled liquids dynamics |l5f| . Then, 
we investigate the statistical properties of the components of VF(x(£)) in this special frame. We find that these 
components are distributed according to a self-similar form, determined by the value of the exponent k of the energy 
decay. 

We develop up to the second order the potential around the actual position of the particle x(t). 
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Q(x) = V(x(t)) + ^(x 4 - Xl (t)) • aiV(x(t)) 
+l/2^%y(x(t)) • -**(<)) • (x,-- Xi (i)), (5.1) 

Wc define an orthonormal frame of eigendirections {e^i)} in which the hessian TLij{t) — didjV(x(t)) is diagonal. Xi 
belongs to the -time independent- interval [0, 2A], so that the corresponding eigenvalue of Q is just Xi — S(t). 

We follow "adiabatically" the eigenvectors {ex^t)} as the particle moves. A mild assumption is that the {e.^} 
evolve smoothly, provided the levels A^ are allowed to freely cross each other. This choice implies that any ordering 
of the Xi lasts only for a short period of time. The {e,^} define a comoving frame, in which the gradient VV, or 
equivalently the velocity, can be projected. 

-VV( X (t))=J2l l (t)-e x M, (5-2) 

i 

= *(*)• 

There are reasons to consider that the components fi(t) are randomly and evenly distributed, even in the deterministic 
zero temperature limit. First, this randomness reflects the average over the "white" initial conditions. Then, as the 
correlator d 2 jV(x)d 2 jV(x') is exponentially short range correlated, one can suppose that the comoving frame is rotating 
on itself in a chaotic manner, as it does in the spherical (p > 3)-spin model [Q. So, during the particle's motion, each 
component spreads continuously over the N — 1 others directions. 

The sign of 7i(i) itself is irrelevant, because of the arbitrary definition of the frame, invariant under the reflections 
e Xi <-» — e Xi . We claim that 7?(t) has to be preferred to ji(t). On physical grounds, we propose to consider only 
the smoothed quantity |~7i(i)_|, obtained by averaging locally over the few v^/V indices j such that A,; — N~ x / 2 < 
Xj < Xi + N 1 / 2 . This is possible because the mean interspacing between the A^ is 0(N^ 1 ). As N goes to oo, one 
expects py? (t) \ to become a smooth function of Xi, varying only on the scale <5A ~ 1 (although, rigorously, the scale 
of variation is 8X ~ N~ 1 / 2 ), making the dependence in the index i irrelevant. 

The function: 

g(X h t)= [ 7l 2 Wj, (5.3) 

is the distribution of the gradient's coordinates (or equivalently of the instantaneous velocity coordinates) and is a 
central object in the present study. 

In this continuous limit, the two first time derivatives of £ can be expressed with the help of the density p(X) and 
the distribution g(X,t) as: 

^)=-Ew))-W)). 



dX p(X) g(X,t). (5.4) 

£(*) = E • W x (*)) ■ W x (*))> 



dXp(X) (X-S(t)) g(X,t). (5.5) 
In |7) was already noticed that, due to the algebraic decay of the energy £ (t) — —2 + 1.08 • t~ K , the ratio £(t)/£(t) 



1/t. From section HI, we know that S(t) ~ t 0,67 , which implies £(t) -C S(t) ■ £(t), and thus: 

f dA p(X) X ■ g(X, t) = £{t) + S(t) J dX p(X)g(X, t), (5.6) 
S(t) fdXp(X)g(X,t). (5.7) 



The first moment of g(X,t) p(X), is proportional to S(t), suggesting a self-similar scaling form for g(X,t), valid for 
t ->■ oo and T = (Figure §: 

g(X,t)=T(t)G(J^-Y (5.8) 
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The knowledge of the other moments of g would be useful to confirm equation (5.8), but unfortunately, they are very 
difficult to compute, and are no more given by the next derivatives of £ . 

As t increases, only the smaller Ai keep any relevance, and the densi ty p is well app roximated by its A ~ equivalent 



tt 1 (2/A) 3 / 2 -\/A- In this limit, the loss of energy rate becomes, from (5.4) and (5.S) 



S(t)=-^-£(t)cx-T(t)-S(t) 3 / 2 . (5.9) 



The knowledge of the exponent k of S(t) ~ t K (section III) fixes the prefactor T up to a constant, to: 



r =5 (2-«)/2«_ ( 51() ) 

Our section || suggests K is very close to 2/3, which would imply r cx S. 

The next momentum of g(A, t)p(X) provides information on the time correlations of the unit vector w(t) of the 
particle's trajectory. One the one hand, 



(3 t W(x(t))) 2 = J2 E -d 3 V{x{t))d l3 V{x{t)) , 



J2-dkV(x(t))d lk V(x(t)) 




= d i v ■ d a v ■ d *kV ■ d k v, 

dXp{X) (X~S(t)) 2 g(X,t). (5.11) 

With the scaling form for g(X,t), the right hand side is of order T(t)S(t) 7 ^ 2 . On the other hand, we perform a 
decomposition — W(x(t)) = M(t) w(i). The norm M{t) equals (— f(t)) 1 / 2 , and w(i) is the unit vector, tangent to 
the trajectory. The following equality holds: 

(d t VVf = (d t M) 2 + M 2 ■ \\d t w\\ 2 . (5.12) 

This sum is clearly dominated by M 2 ■ ||9tw|| 2 , with M 2 ~ T ■ S 3 ^ 2 . The unitary vector rotates, regardless to the 
actual value of T(t), at a rate ||<9tw|| ~ S(t). One expects the "director" w(t) to have changed its orientation after 
a typical time <S(£) -1 , which looks like a "persistence time" for the trajectory of the particle. Consequently, the 
motion of x(t) crosses over from a "ballistic" regime ||x(£ + St) — x(t)|| ~ M 2 ■ (5t) 2 ; St -c S^ 1 to a diffusive re gime 



x(t + St) — x(t)|| ~ D ■ St; St S 1 . The existence of a diffusive regime is here inferred by the expansion ( 2.21 



2.23), valid only for exponentially correlated disorder, and by no means generic. 



VI. A SHORT TIME, QUASI-STATIC APPROXIMATION 

We investigate here the breakdown of the fluctuation dissipation theorem, in the zero temperature limit. The 



fluctuation-dissipation violation is measured by the function T{t,t') (cf 2.17). We propose here a model for the short 
time evolution of T(t,t'), and show that its predictions are in good agreement wit h th e findings of section [HI 



We approximate locally the potential V(x) around x(t) by a quadratic function (5.1), which may can be considered 
as constant provided we restrict ourselves to a time separation t — t' small enough. One can always find a coordinate 
system {y^} such that this quadratic potential reads: 

Q(y) = Q(0) + l/2^(A 4 -5)- 2 / 2 , (6.1) 

i 

where the coordinates y must not be confused with the original coordinates x of the rclaxational motion. 

This section aims at demonstrating that when a particle diffuses, or relaxes in such a parabolic potential, then a 
characteristic time tf scaling like <S -1 arises, which turns out to be the time scale along which the function T(t,t') 
departs from its equilibrium value 0, i.e the fluctuation-dissipation violation characteristic time. 
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We consider a particle moving on the potential (6.1), starting at to, and define the time difference r = t — to- The 
intermediate steps of the calculation make use of r, to, while the final results are expressed in term of t, t' in relation 
with the original out-of-equilibrium relaxation. 

Let us consider the same local average [yf (r)J as in equation (5.3). The are related to the gradient's coordinates 
by J i (r)=-(X i -S)-y i (T): 

r 2/ \ I ff(A,t + r) 

M T )J = ( A _5)2 ■ ( 6 - 2 ) 

The initial conditions [~yf(T = 0) | are given by g(X,to)- One computes the fluctuation dissipation violation T(t,t'), 
when the quadratic potential ( |6.1| ) does not evolve with time (<S fixed once for all), and with initial conditions arising 
from a realistic distribution g(X, to) = T ■ G(X/S). 

Vt{r) = W (0) e-^- 5 )-W. (6.3) 
The distribution <?(A, t + r) evolves like [y? (r)J • (A — S) 2 . One has, far all t > t : 

g(X,t +r)=g(X,t ) e - 2 < A - 5 )'M; 

9 t5 (A,t) = -2 g (A,t)-(A-5). (6.4) 

The usual response r(t, t') = iV _1 ^\ Syi(t)/SQ(t'), and correlation 6(t, t') = A^ 1 (y» (i ) — yi(t')) 2 functions reex- 
press in terms of g(X, t'): 

r(t,t') = j dA p(A) e -( A --s) (t-t') ; ( 6 .5) 

b(t, = / dA p(X)g(X, t') ■ ; (6.6) 



d t ,b{t,t') = -2 / dA p(A) 5 (A,t') 



1 _ e -(A-S)-(t-i')\ 
^ j ; 

= -2T(t,t'). (6.7) 
By inserting t7(A, t' = to) = rG(A/<S) in ( |6.7[ ), one deduces the short time, r <C S^ 1 , value of T(t, t'), 

T(t,t') = (t-t') -r-5 3/2 , (6.8) 

and the intermediate time r ~ one, 

r(t,t')=5 3 / 2 0>o(5-(t-t')), 
6(t,t') = -2r-5 1 / 2 -$ 1 (5-(t-t')), 

T(t,*') = ~(5 -(t-O). (6-9) 



^o^i are scaling function presented in appendix Equation ( |6.9| ) shows that 5" 1 plays the role of a characteristic 
time for the onset of the effective temperature T. 

Assuming now the very likely value k = 2/3 and r = S, one finds that T becomes an order one quantity after a 
time tf ~ One concludes that the characteristic time scale t/ should scale like t/ = t Q oc t K , and a = k = 2/3. 

Our numerics (Figure |]) lead to an estimated value a ~ 0.64, while k ~ 0.67. While we haven't proved that k is 
actually 2/3, we find the agreement satisfactory, and believe that the above picture describes correctly the first stage 
of the breaking of the "fluctuation-dissipation r elat ion at zero temperature" . 

How long can the quadratic approximation (5.1) accurately describe the original relaxational process ? As S(t) 
decreases algebraically, the necessary time St to have |<S(f + St) — S(t)\ ~ S(t) is t itself. More seriously, we have seen 
in the previous section, that the unit vector of the trajectory x(t) changes with the time scale <S -1 ~ t 2 / 3 . As this 
change is somewhat related to the frame's chaotic motion, we deduce than S^ 1 must be an upper limit of validity 
of the quasi-static approximation. Finally, the relaxation on the saddle becomes ill-defined when S ■ (t — t') 1, 
due to the exponential divergence of the functions r(t,t') and b(t,t'), given by the equations (3.5, [6^ ). We arrive 
to the conclusion that this quasi-static picture is not valid beyond times much greater than S~ , but provides a 
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strong presumption in favour of tf(t) ~ ~ i 2 / 3 , in good agreement with our numerical findings (section [II and 

Figure |j) . 

Let us close this section by computing the typical distance covered during a time interval t — t' < with a 

gradient coordinates distribution g(X, t') = T ■ G(X/S). 



b(t,t') = rs 3/2 • (t-t') 2 

= S 5/2 ■ {t-t'f if k= 2/3. 



(6.10) 



For a time interval {t-t') ~ S' 1 , b(t,f) - r-5" 1/2 < 1, becoming 6(t, t') ~ 5 1/2 for k = 2/3. As 5 1/2 tends to zero, 
the characteristic time t^ of the evolution of b(t,t') in the aging regime, is necessarily much greater than tf ~ <S _1 . 



VII. A DYNAMICS RESTRICTED TO THE DOWNHILL DIRECTIONS 



This section shows how the above approach describes the long time aging regime. 

The equation ( |5.9D has a simple physical interpretation. The only non-vanishingly small components 7$ of — W 
are those corresponding to A^ < S. Only a number N J Q dAp(A) ~ N -S 3 ^ 2 directions i are contributing to — (VT^) 2 = 
— V\ 7? ■ Each one of these 7, has a magnitude of order f 2 ~ I\ As a result, AT • £ scales like — f 2 — N ■ S 3 / 2 ■ T. 



£^-(£(t)-£ d ) 1+1 /-. 



(7.1) 



The relaxation dynamics looks like if it was controlled by the difference £(t) — 

The linear response regime. A constant force T is now applied, uncorrelated to the potential V. Each one of 
its (comoving) coordinate /, is random, time-dependent as the frame rotates during the particle's motion, and has a 
magnitude fi ~ F. We suppose F weak enough to be considered as a perturbation around the relaxational dynamics 
described in section VI. 

At any time, there are "open" , or downhill directions, with Xi < S and "close" , or uphill directions with Xi > S. 
The close directions behave as confining harmonic potentials which prevent the (weak) force fi to drive the particle 
along this direction. The open directions are the one along with the external force drives efficiently the particle 
away. As the particle moves, "open" and "close" directions exchange their role, but the proportion of open directions 
remains proportional to S 3 ^ 2 . 

The force J- induces a displacement x whose components are ~ fi along an open direction and x, ~ along a 
close direction . The average velocity x • .F/||.F|| is given by (6 Heaviside function): 

= N-S 3/2 -F 2 /(VnF); 
= VN F S 3/2 ; 



(7.2) 



that we identify to v N u(t). As a result, one finds a velocity proportional to the number of downhill directions: 

u(xF-S 3/2 . (7.3) 

Inserting the likely value k — 2/3, one finally gets a displacement u(t) — u(t') oc F(\n(t) — ln(f')), well confirmed by 
the numerics (section III and Figure [|). This pure relaxational motion is driven by the components 7^ of — W along 
the open directions, while the external force acts with /, along the same open directions. One expects the linear 
response to hold if f 2 <§C 7 2 but to break down when f 2 ~ 7 2 . This leads to a predicted cross-over time scaling like 
T(t F ) = F 2 , or t F ~ F 4 /^- 1 ^, to be investigated in a forthcoming publication 127 



The diffusive regime. The asymptotic behaviour predicted for b(t,t') is, from equations (2.21-2.25): 



b(t,t') = 



t-t' 



O 



t-t 1 



(7.4) 



One recognises a simple diffusive behaviour, with effective diffusivity t b . From section |IV|, we know that the short- 
time motion (t — t') < S^ 1 ~ tf is ballistic, and that the particles covers a distance of T ■ S~ x / 2 On the other hand, 



our results from section III show that the direction w of the trajectory x(£) uncorrelates itself after this same time 
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S 1 . Using a well-known result on correl ated random walks, and assuming a free diffusive behaviour at intermediate 
times t — t' ~ tb, as inferred by equation (|7.4|), one finds: 



This corresponds to ballistic steps of length (T/VS) ( |6.10 ), and a cross-over time from ballistic to diffusive regime 
equal to As tf ~ <S _1 , 

&(t,i')^(i-t')-r-5 1/2 , (7.6) 

leading to the identification: 

t b oc r- 1 ^ 1 / 2 ; 

5 1/K (7.7) 



oc 



If k is taken to be 2/3, one gets % ~ t' ~ t. 

Let us discuss different reasons to be confident in the scaling tf ~ tb ~ 6>~ 3 / 2 and 5 ~ t~ 2 l 3 . First, a matching 

argument similar to [ p~6|J]~l] ] predicts tb ~ i^ 2 . Then, the result t& ~ £' is in agreement with the conjecture h{t) ~ i 5 . 
This entails a logarithmic growth of b(t,t'), t' fixed, and we have asymptotically (i.e. t,t' ^> 1, and t/t' ~ 1) a free 
brownian motion in logarithmic time. 

b(t,t')=5-(\nt-\nt'). (7.8) 

This makes exp(—b(t,t')) as well as r(t,t') decaying as a power law. While we have no demonstration of that, we 
think that a p ower -law decay of the memory function f(b(t, t'))r(t, t') is necessary for the "fine tuned" aging solution 
of the system (|2.8|). Asking for a power-law decay f(b(t,t')) in turn fixes k to 2/3. 

Finally, if K = 2/3, tb = <S~ 3 / 2 , and the characteristic times for the linear response regime u(t) ~ F/tb, and for the 
diffusion regime b(t,t') ~ (t — t')/tb are the same, which is consistent with the persistence of an "Einstein relation" 
at the beginning of the aging regime. 



VIII. CONCLUSION 



We have proposed a geometrical description of the mean-field relaxational dynamics of a particle, for a subclass of 
short-range correlated disorders. We have restricted ourselves to the isolated case, and to the driven case in the linear 
response regime. 

A numerical integration of the mean-field equations gives evidence of a power-law decay of the dynamical energy 
with an exponent k numerically close to 2/3. We also found evidence of a logarithmic growth bit, t') ~ Inf consistent 
with the conjecture hit) ~ t s for the reparametrization function h. 

The exponential correlator makes it possible to compute the density of eigenvalues of the hessian Ti associated to the 
random potential, and we were able to predict the correct value {i.e. —2) of the dynamical energy Introducing a 
comoving frame, reminiscent from the INM frame of a supercooled liquid, we derive an expression for the distribution 
g(X, t) of the components of W(x(i)) . This expression is g = T G(X/S), where —S{t) is the (time dependent) lowest 
eigenvalue of Spec(7Y). 

For reasons exposed in section |Vl4 namely the consistence with h(t) ~ t 5 , the requirement that /(&) is likely to 
decrease as a power law, and acknowledging the numerical estimate of k, we believe that k is indeed equal to 2/3. 
This leads to the following predictions: 

(1) F oc S, and a typical gradient coordinate is, along a downhill direction, |tJ ~ ~ £ -1 / 3 . 

(2) From a short time, harmonic expansion of the particle's motion (section |V|) the characteristic time tf leading 
to the appearance of an effective temperature goes like tf ~ S^ 1 ~ t 2 / 3 . 

(3) The characteristic time at the beginning of the aging regime is tb ~ S~ 3 / 2 ~ t. Both linear response u ~ F/tb 
and diffusion b(t,t') ~ (t — t')/t b are controlled by it. 

We conclude that the aging mechanism of this model comes from a simultaneous decrease of the number of downhill 
directions (going like N ■ S 3 ^ 2 ~ Nt^ 1 ) and of the typical gradient component |7i| ~ t -1 / 3 . 

We predict also that the effect of a constant force brings about a dramatic change in the dynamics after a time 
tp ~ F~ 3 , reaching a out-of-equilibrium but stationary regime p7[ . 
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APPENDIX A: THE MSR ACTION 

The action leading to equations (|2.8| ) is: 

,00 ( n N ) 

S[x,ix] = / dt \ -T V(^) 2 (i) • - ^ 1/2 • + V i% (i) • % (*) \ 

J ° { 3 = 1 3=2 J 

roc ( N N } 

+ dtds {f\d(t,s))J2ix j (t)-ix j (s) + 4f''(d(t,s))r(t,s) V) *%(*)■ (&i(t) - > , (Al) 

J ° { 3=1 3=1 J 

and the expectation value of an observable 0[x(t), ix(t)], averaged over the disorder, is given by: 



(0)= Vx[t] Vx[t] O exp(-S). (A2) 



APPENDIX B: THE SPECTRUM OF THE HESSIAN H 



We consider 6(V(r) — Vq) Spec(7i)(r) for any arbitrary potential. Vq = NE is fixed, and th e Tij j = dfjV(r) are 
N(N + l)/2 gaussian random variables. The correlations among the Tlij are listed in equation ( |4.6| ). We define the 
self-averaging quantity V = N~ x J2i duV(r) so that: 



N(d vl V(r) - vy = 8/"(0) - 8/"(0)/iV, 



N(d u V(r) - V)V(r) = 0, 



N(d u V(r) - V){d n V{v) -V) = -8/"(0)/iV, 

= + O(A^- 1 ). (Bl) 

The hessian is now Tiij = T>6ij + 0.'^. W is a matrix of independent gaussian centred random numbers. The diagonal 
elements are slightly correlated (of order 1/N 2 ) and have a different variance than the off-diagonal elements, but this 
does not prevent the Wigner result to apply and the spectrum of 7i' is a centred semi-circle of width A = 4-%/ /"(0). 
The determination of T> follows from the fact that £ and T> are gaussian distributed, with correlations: 

Ne s = f(0); 

NV^ = 4(^±- r ) f(0); 



NV ■ E = 2/'(0). (B2) 
The joint probability distribution of E and T> is: 



£ 2 (V - a£f 



/(0) 



4 (/(o)/"(o) - /'(o) 2 ) + 4/"(0); 



/(0) W V . , KJ , N 

2.f(0) 
/(0) 



(B3) 



Fluctuations of 2? are of order A 1 / 2 around its saddle point value 2/'(0)//(0) x It follows that Spec(7Y) is a 
semi-circle of radius A shifted by an amount V = 2/'(0)//(0) x £. 
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Let us consider now 



S(V(r) — Vq) 5(V(ri) — Vx) Spec(W)(r). This simple average measures the non-locality of Spec(H), i.e. its dependence 
on the values taken by the random potential V(r') around r. 

The rotational invariance of the above average is broken, and r — r' plays a special role. We relabel hereafter the 



direction 1 to coincide with r — r', and define b = ||r — r'|| 2 /7V. Correlations are now, in addition to (4.6) 



d ll V(T)V{v') = 2f'{b) if z > 2; 
dnV^Vir') = 2f'(b) + 4b f "(b); 
d l] V(v)V(r') = 0. 



(B4) 



V is again defined as iV" 1 i duV(r) and H'^ ^ >2 is equivalent to the above situation: independent, centred, 
gaussian random components, and the spectrum is a centred semi-circle. Adding one row and one column of random 
independents elements to H'^ ^ >2 must not change the density profile of eigenvalues. This is because this eigenvalue 
distribution is a fixed point under the change N — > N+l, as argued in the cavity approach of the problem. A possible 
trouble come from the single component d\\ V — T> which does not average to 0, but this does not alter the final result 
more than by a single isolated eigenvalue. 

It is possible to show, with the help of a formal field theoretical approach, that the correlations ( |4.6| ),(B4) indeed 
lead to the ordinary N — > oo saddle point for Spec(W), i.e. a semi-circle law of radius A. 

The computation of T> follows closely the lines of the previous paragraph. We found that if £ = V(r)/N, £' = 
V{r')/N and b = ||r- r'\\ 2 /N, then: 



V 



f(0) + /'(&) 
/(0) + /(6) 
f(0) - f'(b) 
/(0) - f(b) 



{£ + £') 
■{£-£'). 



(B5) 



For a generic correlator, there is an explicit dependence on £' ("non locality") while for an exponential correlator 
/ = cxp(— y), the above formula reduces to V = 2/'(0)//(0}x £. 

This suggests that the determination of Spec(7i) from ( f4.2| ) is a complex problem and the simple behaviour (4.4) 
fails for a generic /. 

The exponential correlator, however, has a strong property. The average: 



N 



/"(0). 



(B6) 



is at the order N 1 . This means that, while b is strictly positive, there is no possible fluctuations of T> around 
2/'(0)//"(0) x £ . 



Repeating the argument for (4.9), n fi nite , shows that Spec (7i) depends only on V(r) and not on its local environ- 
ment. Thus, we argue that the average (4^2) is given by (O, 4^), as announced. 



APPENDIX C: THE QUASI-STATIC PICTURE 

From equation (|6.5|), we derive the expression for r(t' + r,t'). 

e ST Ii(Ar) er Ar 



r(t'+T,t') = 2- 



At 



(CI) 



where Ii is the first kind modified Bessel function. The short time expansion of (6.7) is: 

d t ,b(t' + r, t') = -2r / dA p(X)g(\, t'); 



= -2r(£(t')) ; 
= -2r TS 3/2 . 



(C2) 



16 



In the intermediate time separation regime, r is of order S 1 . The integral is dominated by A <~ S and cut off by 
g(X,t') for XS 3> 1. p(X) can be replaced by its A — > equivalent. 



r(t' + r, t') ~ v^e^^Ar)- 3 / 2 ; 

= S^y^A" 3 / 2 ) • (Sr)- 3 / 2 e 5T ; 
= 5 3 / 2 $o(5r). 

d t >b{t' + t, t') = -2TS 1 ' 2 / du (2/A) 3 / 2 tt- 1 Vu 

I _ „-St(«-1) 1 

xG(«) ' 



The effective temperature behaves as: 



reducing to 



u- 1 
-2r5 1/2 $i(5r). 



T(<+r ' t)= UJ^)' 



m0^(5'(t-0), 



(C3) 



(C4) 



(C5) 



(C6) 



if k = 2/3 andT = 5. 
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Captions 



FIGURE 1. Parametric plot of the integrated response b(t, t') vs 7Z(t) — J* ds r{t, s) at zero temperature, for time 
steps h — 0.025, 0.5, 1.0, 2.0. The horizontal part corresponds to the short time regime, with T — ► 0. Then, the aging 
regime is the straight line with a slope X^ 1 ~ 0.21, to be compared with the theoretical value 1/Xqfd = 2. The 
inset shows the derivative X _1 (6) = db/dlZ(t), stepping from to 2.1. 

FIGURE 2. Dynamical energy £(t) + 2 vs time, in log-log coordinate, for h = 0.025 and h = 0.2. The power-law 
decay is unambiguous, and a fit to k — 0.67 has been done between the two vertical arrows. 

FIGURE 3. Logarithmic derivative — k of the energy £(t) + 2, for h — 0.025,0.5,1.0,2.0. The curve is noisy as 
£(t) — » —2. The straight line stands for k = 2/3 which we believe to be its exact value. Curves for h — 0.025,0.05 
seems to tend to 2/3 from above, while h = 0.2 seems to tend to 2/3 from below, k ~ 2/3 is well realised for h = 1. 

FI GUR E 4. The characteristic times tf(t) oc t a and t a (t), a = 0.55 and a — 0.45 determined from equations 
( |3.2| ) and (3.3). Results are shown for the time steps h = 0.1 and h — 0.2, and the finiteness of h is visible at small t. 
A numerical estimate of a is 0.64 between the first and last vertical arrows. The exponent of t Q =o.55 is close to 0.93 
while we expect 1, and t a= o.45 should saturate to a constant. 

FIGURE 5. The functions exp[— b(t, 0)], exp[— b(t, t—tf(t))] vs t, for h — 0.2 and h = 0.1, in logarithmic coordinates. 
The functions exp[— b(t, 20)] vs t — 20 and exp[— b(t, 40)] vs t — 40. Here, 20 and 40 are waiting times. The behaviour 
of b(t, 0) and b(t,t — t /(t)) is doubtless logarithmic. The slopes of exp[— b(t, 0)], exp[— b(t, t — tf(t))], on this figure are 
respectively —1.10 and —0.42. According to the predictions of pOj , exp(— b(t, t')) tends to h(t) — h{t') for t,t' » 1; 
t/t' finite. The curves exp[— b(t, 20)] and exp[— b(t, 40)] tend to imitate exp[— b(t, 0)], with a delay. 

FIGURE 6. A test of the linear response of the displacement u(t). We plot u(t)/F, as a function of ln(i), for 
F = 0.05, 0.1, 0.2 and 0.4; h = 0.2. As F — > 0, the curves are indistinguishable from the integrated response 7£(i). A 
departure from the straight line signals the breakdown of the linear response, as the particle acquires a finite velocity, 
dependent (non-linearly) on the force. The suggested behaviour of u(t) is thus: u(t) — F ■ (c^ + C4 hi(i)). 

FIGURE 7. The density p(X) of eigenvalues A — S (on top). A sketch of the self-similar distribution g = SG(X/S), 
assuming r = S (bottom). The tail of G goes to as X/S — * 00 , in the asymptotic limit 5^0. 
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